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Abstract: We show that boundary contributions of BCFW recursions can be interpreted 
as the form factors of some composite operators which we call ’boundary operators’. The 
boundary operators can be extracted from the operator product expansion of deformed 
fields. We also present an algorithm to compute the boundary operators using path integral. 
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1 Introduction 

It is well known that BCFW recursion relations [1, 2] (see also reviews [3-5] ) fails to 
capture a piece of scattering amplitude, usually called boundary contribution, when the 
scattering amplitude does not have the desired vanishing large z-behavior under a given 
BCFW deformation. After years of development, there are several methods proposed to 
solve this problem, such as introducing auxiliary fields to eliminate boundary terms [6, 7], 
analyzing Feynman diagrams to isolate boundary terms [8-10], expressing boundary terms 
as roots of amplitudes [11-13], and collecting factorization limits to interpolate boundary 
terms [14]. However, these methods are applicable only to limited types of theories, and 
the evaluation of boundary contributions in more general theories remains a difficult task. 
Recently two ideas have been proposed to solve boundary contributions. In [15-17], multiple 
steps of BCFW-deformations have been used to track boundary with certain poles one 
by one, while in [18], multiple legs have been used to do the deformation. These two 
methods can be applied to broader theories although the pure polynomial part of boundary 
contributions remains a challenge. 

In [16] it was found that with proper choices of deformations boundary contributions 
satisfy similar recursion relations as scattering amplitudes. Based on the observation, a 
systematical algorithm was proposed in which scattering amplitudes are computed with 
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multiple steps of recursions. The first step is the usual BCFW recursion relation, which 
leaves out a boundary contributions^ The second step is a recursive computation of 

which leaves out a new boundary contribution and so on. A fe-step recursion 

yields the complete amplitude if vanishes. Especially it has been proved that 

scattering amplitudes of 4D renormalizable theories can be computed by (at most) 4 steps 
of recursions. 

The multiple step recursion method could be further refined in several directions. First, 
in [16] the large z-behavior is determined through a careful analysis of Feynman diagrams, 
which was already cumbersome for the examples considered therein. More convenient and 
systematic methods are needed for the discussion of models with more complicated matter 
content and interactions. The next problem is under what condition could we decide the 
large ^;-behavior of a deformation through some direct tests of simple amplitudes. Suppose 
under a certain deformation, the large ^;-behavior of any amplitude with number of legs 
n < N is 0{^), then it is tempting to speculate that the same goes for any n. However, 
this is not always true, at least if the Lagrangian contains a m > N point vertex which 
only contributes to amplitudes with n > m. Last but not least, in many theories even 
after exploiting all possible deformations, the last boundary contributions 
do not vanish. The first way to solve the problem is, of course, developing a supplementary 
method which computes directly. Another avenue to attack this problem is to 

explore whether ij(i"-"--2|«.] ^^n be discarded at some limit(for example suppressed by some 
large energy scale in an effective field theory) even if ^ q, so that the recursion 

terms would be a good approximation of the complete amplitude. 

Moving forward in any of the directions above requires a more rigorous understanding 
and transparent formulation of boundary contributions. In this paper, inspired by operator 
product expansion (OPE), we relate boundary contributions to some composite operators 
which we call ’boundary operators’. Boundary contributions can be expressed as the form 
factors involving the corresponding boundary operators and the undeformed external states. 
The employment of this new formulation helps to solve all the questions above. For the first 
two questions, the large z-behavior can be read directly from the OPE. The vanishing of 
boundary operator is equivalent to the vanishing of boundary contributions for all amplitude 
under the corresponding shift. For the last question, the last boundary contribution can be 
read directly from the corresponding boundary operator, since there is only one undeformed 
external leg left. 

Boundary operators serve as a bridge connecting on shell and off shell quantities to¬ 
gether. During the last two decades the computation of scattering amplitudes has been 
facilitated by the newly developed unitarity based methods, but these new methods usually 
do not apply to off shell quantities like form factors and correlation functions. So although 
not the focus of the current paper, our result will open a new window into the computation 
of form factors and correlation functions, and can be seen as part of the effort devoted to 
extract all physical quantities from scattering amplitudes. 

means the boundary contribution to the amplitude under the (l|n] deformation. means 

the boundary contribution to under the (2|n] deformation. 
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The structure of the paper is as follows. In section 2, we briefly review the bound¬ 
ary contributions, and discuss relations between boundary operators and operator product 
expansion(OPE). In section 3, we present a path integral derivation of the boundary oper¬ 
ators. In section 4, we give some examples to demonstrate the computation of boundary 
operators, and in section 5, we give some concluding remarks. In the appendixes further 
examples are presented. 

2 Boundary Contributions and Operator Product Expansion 

A well known fact about BCFW recursion is that the large 2 ;-behavior of amplitudes under 
a given BCFW-deformation usually only depends on the type of two deformed legs, not 
the type or number of undeformed legs^. For example, the large z-behavior in Yang-Mills 
theory can already be seen from 3-point or 4-point MHV gluon amplitude. This interesting 
’locality’ property of BCFW recursion is also inherited by multiple step recursions in [16], 
and by (T, X) multiple line shift discussed in [18]. 

This ’locality’ property of large z-behavior is most naturally described in the language 
of the operator product expansion (OFF). More explicitly, in a quantum held theory, the 
product of two nearby operators Oj{x) and Oj{y) can be expanded by a basis of local 
operators: 

Oj{x)Oj{y) = Y,Cjj^{x - y)OK{y) (2.1) 

K 

In momentum space, (2.1) is translated to 

Oi{pi)Oj{p2) = Y,Cij^{pi)Ok{pi +P2) (2.2) 

K 

where the condition that x and y are nearby in coordinate spaces is equivalent to the 
condition that pi and p 2 are quite split in momentum space. If we set pi to ki + zq, and 
P 2 to kn — zq, then the OFF holds when z is large: 

Oi{ki + zq)Oj{kn - zq) = '^Cjj^{ki + zq)OK{ki + kn) (2.3) 

K 

If the momentum of an operator is not shifted by ^zq under the BCFW shift, we will 
call it a soft operator. And if the momentum of an operator is shifted, we will call it a 
hard operator. The total momentum of two hard operators with approximately opposite 
momentum is soft, and (2.3) tells us that the effect of these two hard helds must be described 
by some local soft operator. 

The r.h.s. of (2.3) depends on z only through some c-number functions Cjj^{ki + zq). 
Expanding coefficients around z = oo, 

+ zq) = ^ C^jj^z^ , (2.4) 


^Of course, as pointed out in the introduction, sometimes boundary contributions of a shift vanish 
accidentally for amplitudes with small n. 
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then we have the following asymptotic expansion of (2.3) at large z, which we call G{z)^ 


G(z) = Oi{ki + zq)Oj{kn - zq) = '^ Qiz\ 

i 


(2.5) 


in which Qi = C^jj^Oxiki + kn), and each Gi is a local operator. 
In the language of form factors, (2.5) becomes 


Oi{ki + zq)Oj{kn - zq) $(^ 2 ) • • • ^{kn- 2 )^ = X] ^{^ 2 ) ■ ■ ■ ^{kn- 2 )^- (2.6) 


If we set Oj and Oj to be two fundamental fields ^{ki + zq) and ^{kn — zq), after LSZ 
reduction'^, (2.6) becomes an n point scattering amplitude ^(^(fci),-- - ,^{kn)) under the 
(l|n] shift, and the boundary contribution under this (l|n] shift can be written as 


= {0\gomk2)---Hkn-2)) . (2.7) 

In this paper, we are mainly interested in Go^ and for later convenience we will also 
denote this operator as and call it the boundary operator of the BCFW shift. 

Similar steps leads to the boundary operator of a multiple step shift 

= {0\O^^-P\^\Kip)Mkp+i) • • • Hkn-2)) . ( 2 . 8 ) 

where Kip = ki + k 2 + ■ ■ ■ + kp. The only difference is instead of being a fundamental field, 
G^j{kR — zq) in (2.3) is now the boundary operator of the previous shift 


3 Derivation of the Boundary Operators 

In the last section we found that boundary operator is the O(z^) order of the OPE of two 
deformed helds. In this section we present an algorithm to compute the boundary operators 
via path integral. We split each field into a hard field and a soft field, and the OPE can be 
interpreted as the Green’s function of these two hard helds, and can be evaluated exactly 
at tree level. 

The boundary operator of a multiple step BCFW deformation corresponds to the OPE 
of a fundamental held and the boundary operator of the previous deformation. In order 
to have a uniform description of boundary operators of the hrst deformation and the later 
deformations. In subsection 3.3, we give a uniform description of boundary operators of 
the hrst deformation and the later deformations by introducing the boundary helds. The 
boundary operators of later deformations correspond to the OPE of a fundamental held 
and a boundary held. 

^In the next section we will see G{z) can be interpreted as the Green’s function of two hard fields in the 
background of soft fields. 

“^This means strictly speaking Oi{kL) and Oj{kii) are not simply the fundamental field. Instead they 
are fundamental fields with LSZ reduction. 
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3.1 A Path Integral Representation of the Bonndary Operators 

In this subsection we present a path integral presentation of boundary operators. Consider 
the following n-point correlation function in a general quantum field theory: 

Gn{z) = mkimk2)---^{kn)) = I D^exp{iSm^ki)Hk2)---Hkn). (3.1) 

In order for the OPE (2.5) to hold, the momentum of the deformed particles $(A:i) and 
^{kn) must be very large. Following Wilson’s idea, we can split the field into a high energy 
(hard) part^ <I>^, and a low energy (soft) part (which is still denoted by <I>), 

<I> ^ <!> + . (3.2) 


Under the splitting, the action S'[<h] becomes 

5 [$ + = 5 [$] + $] + $] + ••• 


(3.3) 


where $] contains terms linear in and S'^[‘I>^, <I>] contains terms quadratic in 

and so on. Now we choose the energy scale A to satisfy \zq\ ~ A S> \ki\. Shifted lines 
(propagators or external legs) carry momentum p ± zq, thus correspond to the hard field 
while undeformed lines correspond to the soft field $. 

A typical tree diagram containing two hard external legs and multiple soft external 
legs is shown in Figure 1. Obviously all hard fields he in the line connecting the shifted 
momenta ki + zq and kn — zq, and each vertex in the diagram contains either zero or two 
hard fields. This means only S'[‘I>] and S'^[‘I>^, $] contribute to Gn{z) at tree level. 


^\kL + zq) 


\ 1 / 1 

\ 1 / 

\l/ ^, 1 

\ \l/ 
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I / 
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Figure 1: A diagram contributing to the deformed amphtude(correlation function). Solid 
lines represent hard <I>'^ fields, while dashed lines represents soft ‘h field. Each vertex in the 
diagram contains zero or two hard fields. 


Then the shifted correlation function can be evaluated as 


Gniz) = j exp {iS[<^] + $]) • • ■ ^>^-1 

= j Z)$exp(z5[^>]) j T)^'^exp(z5^[^'^,$])$f^> 


(3.4) 




®In [19], fields are split into background and a perturbation parts in order to assess boundary contri¬ 
butions. Eq. (3.2) serves as an interpretation of this splitting, and we are indebted to Mingxing Luo for 
pointing this out. 
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where for simplicity we used the notation 

+ zq), = ^^{kn - zq), = ^ki) . (3.5) 

In (3.4), only terms in the square brackets depend on z, which is the Green’s function of 
two hard helds in the background of soft fields, and will be denoted by 

G{z) = -ij exp CD]) ( 3 . 6 ) 

and can be interpreted as the OPE of and After doing LSZ reduction to both 
deformed fields, G(z} will become G{z) in (2.5). 


3.2 Evaluation of the Path Integral 

In (3.6), <k] only contain terms quadratic in <I>^, and G{z) can be evaluated exactly. 

Now we must be more careful about the components of d>. In 4-D theories, there are real 
fields (gauge helds, Majorana spinors and real scalars) and complex helds (Weyl fermions 
and complex scalars). Suppose the theory has M real helds and N complex helds to 
simplify our notation we combine them into a single held and combine the hard helds 

into a single hard held 


(A 



(3.7) 


The Hermitian conjugation of <I>" is (jy^ (pj ! which can be related to as 

with the matrix 


rpOLj^ 


^ Im 0 0 \ 

0 0 /at 

Y 0 In 0 J 


(3.8) 


as a metric to move indices up. With these notations, the quadratic term in the Lagrangian 
can be written as 
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V^a = 




L, 


(3.9) 


in which is a Hermitian operator. Following the standard procedure of computing 

generating functions, we add a source held J“ for 77", 

= J + (3.10) 

The source held also satisfy J" = J^j^. It is easy to see that in (3.10), jXh°^ = 
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Shifting —>■ — {T> and integrate over the hard field, we find 

Z^[J, j DH exp {iS^[J, H, $]) 

=Z'^[$]exp j d^xd^yJlix) {V-^y (3.11) 

=Z'^[$]exp j d^xd^yJlix) (x, y; $) jj(y)^ , 


in which Z^[‘l>] = detD(<l>) 
index fd using the T-matrix: 


, and in the last line we have switched the positions of the 


Notice T^js can be interpreted as a metric, and we will consider 'D~^ (as well as P, 'Dq, l/)as 
a tensor whose indices can be lowered or raised by or 
From (3.6) we have 


G^f^{z) = - z / PPexp $]) H-{x)Hf^{y) = 

J dJa{x)dJp{y) J-O 

=Z^[$] {V-y^ (x,y;$) 


(3.13) 



2(a) 



Figure 2: The first two diagrams contributing to Z^[$]. 


Now we analyze these two factors in (3.13). In general, the operator P(<h) can be 
decomposed into a free part® Dq and an interaction part V, 

P%(<l>) = (Po)% + 0(‘^)- (3-14) 

The term in (3.13) is a loop effect and can be set to one since in this paper we mainly 

discuss tree amplitude. To see this, first take the limit A ^ 0 and can be expanded 


as 


^ detP($) 
det Dq 


1 - i Tt{D^W) + i Tt{D^^VD^^V) + ^Tr2(p-V) + ■ • • 


(3.15) 


®For example the free part Do is — md for a scalar field, and i(^ + m) for a Dirac field. 
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The first term 


det D, 


0 in (3.15) is a constant factor independent of external fields <h. 

The second term in the parenthesis, —^Ti{Dq^V) corresponds to Figure 2(a). The third 
term, Dq^V), corresponds to Figure 2(b). And the fourth term corresponds to 

a Feynman diagram with two copies of Figure 2(a). With non-zero A, the path integral in 
Z^[^\ cannot be evaluated exactly, but the analysis above suggests 




1 - ^TvAiDr^ 


1 . 


0 ^) + 


-1 


1 


VD^^V) + -TrliD^^V) + 


(3.16) 


in which TrA means the loop integration is over the hard region |/| > A. In all, the Z^[<h] 
term is a loop effect which do not contribute in a tree level calculation, and from now on 
we will set it to one. 

Now let us turn to the factor (P“^)“^ Using the decomposition (3.14), the 

operator Vap can be formally inverted as 


P-H4>) = {Do + U(4>))-i = {Do[l + D^^V{^)]) ^ = [1 + 

oo oo 

k=0 k=0 


Now we move to the OPE G{z), for which we need to do LSZ reductions for fields 
H{ki + zq) and H{kn — zq), i.e., we need to multiply the inverse of propagators and mul¬ 
tiply corresponding external wave functions ei,en. This action is given by the operator 
e^^.{—iDo)°^''l^H^\ Putting it back to (3.17) we get the expansion^ 


g(z) {P-'(-i?) - Oo-')'’*'’- 

= - ‘i/L [Oo - »o"') Oo]“‘“" 

oo 

fc =0 


(3.18) 


From (3.18) we see that it is most convenient to derive Q{z) from T) with two upper indices, 


^ ^ 6^L (3.19) 

It would be worthwhile to stress again that G{z) is the complete OPE which contains all 
z-dependence, not only the 0{z^) order. The 0{z'^) order of G{z) is the operator which cor¬ 
responds to the 0{z'^) order of the deformed amplitude. For example, the {g~{ki)\g^(kn)] 
deformation in Yang-Mills theory is 0{z^), and the 0{z^) order of the deformed ampli¬ 
tude equals to the form factor of the operator Gs- More details on the large z behavior of 
Yang-Mills amplitudes can be found in Appendix B. 

^We subtract a from in (3.18) because after canceling one Do, another Do acts on on-shell 

wave function ei, e„ will give zero. 










3.3 Boundary Operator of the Second Deformation 

In the previous subsection we discussed the derivation of boundary operators of the first 
deformation (l|n], and in this subsection we will discuss the boundary operators of the 
successive deformations. The boundary operator of the fc-th shift (12 • • • fc|n] is given by 
the OPE of Hk and the boundary operator of the previous shift (p(i 2 -fe-i|n] example, 
when k = 2, 

G(z) = -i j DHexp {iS^[H,<l>]) (3.20) 

The in (3.20) is a composite operator, and in order to evaluate the path integral, 

we will introduce a boundary fields* which allows us to rewrite (3.20) as a correlation 
function of two fundamental fields. The boundary fields T couples to the (Hermitian 
conjugate of) boundary operator of the previous deformation, 

L2{H, T, $) = L($) + + o(ih]jr. (3,21) 


Then G{z) can be rewritten as 

G(z) = -i j DHDTexp $]) H 2 T, (3.22) 

where 5^[Pf, $] is defined through the Lagrangian L 2 {H,T,^) by repeating the steps 

described in the last subsection. 

However, unlike the hrst shift, the boundary field is non-dynamical. In order to use 
the result of the last subsection, we will introduce a large mass M to the boundary field 
and assume it is dynamical® , 

L~ = L{^) + jrc)(bn] ^ (3,23) 


and compute the OPE of Pffc+i and following the algorithm of the last section, and set 
Af ^ 00 in the end. 

In practice, since M is large, the term is negligible, then we can equivalently use 
the Lagrangian 

Lm = L{^) + (3.24) 

where M is dehned via Af^ = kM'^. 

Using (3.19), we hnd 


pa/3 


/ Vu 

5<I.t 

50(iN 


s<s>i 

0 -Af2 

-Af2 0 


(3.25) 


®The boundary field J- is an auxiliary field which is introduced so that we could have a uniform description 
of the first and the successive shifts. It will not appear in the final results (3.31). 

®The constant k is needed to offset the mass dimension of J^. 
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where Pff = 1)“^ + and = of + . Again V 

il ' ' f/FsT r^T r^T r^T U ' C> 




be split into a free part T>q and an interaction part V, 

(Dq 0 0 \ 

0 0 -M^ 

0 -M^ 0 y 




fVn - Do ^ 


can 


Vo = 


, v = 


~s^ 

V 


<54-L 

0 0 

0 0 


(3.26) 




Since we are only interested in the shift (<1>“|J^], we only need to worry about the (1,2) 
entry of the matrix V(1 + Dq^V)~^ . Using the large limit, 


p-i _ 

L’q — 


( U>o ^ 0 0^ 


V 


0 0 0 

0 0 0 , 


(3.27) 


it is straightforward to hnd 

[V(l + Po”'V)-'] 12 = (l + (^11 - ^o)^o"') + O(^) (3.28) 

with iy“ = Eq. (3.28) can be further simplihed using the fact that there is no soft 

or .F fields as external fields, so we can drop terms depend on F or .F and finally reach 


[v(i + Vv)-i^ 


112 


1 + UF>, 


-1 


-1 


W 


(3.29) 


Last, we need to multiply this quantity by the wavefunction of H{k 2 ), 

— 1-1 02 




0:2 


l + VD, 


-1 


W^. 


/9 


(3.30) 


Similarly, the boundary operators of the k-th. shift, jg given by the 0{z^) 

order of which has the expression 




Oik 


l + VD, 


-1 


-1 


Ok 

0 5^1 


(3.31) 


4 Examples 

Now we demonstrate the computation of boundary operators by by several examples. We 
will start with the non-linear sigma model, which is the low energy effective field theory 
of pions, focusing on the derivation of T> operator in the presence of multiple fundamental 
fields. Next, we propose a different decomposition of the operator T), which simplifies the 
calculation especially when V{z) ~ 0{z). Finally we discuss the boundary operators in a 
theory with Yukawa and quartic scalar interactions. 
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4.1 The Non-linear Sigma Model 

The low energy effective field theory of pions is described by the following SU (2) non-linear 
sigma, 


L = - 


1 


(4.1) 


2(l + f^)2’ 

where vr = (tt^, vr^, vr^), and T is a parameter with mass dimension -|-1. 

In order to hnd P first we compute the variation of the Lagrangian with respect to 7r“, 


where we have defined 


8L_ 27r“(d^7r)2 

(fvr® ^ F'^oj^ 


TT 




(4.2) 


(4.3) 


In this subsection, we use the convention that if a is inside a parenthesis, then this 5^ 
only act on fields on its right inside the parenthesis. A in a square bracket acts on all 
helds on its right, even outside of the square bracket. 

There are two more subtleties in deriving (4.2). First, should be interpreted as a 
column vector 

/ SL 



5L 

6tt°- 


Second, the operator acts from the left of L, i.e. 

5L 


(4.4) 


5L = 


„ -J-U 


(4.5) 


The next step is compute the variation of with respect to tt°, but this time acts 

5 

67r° 

2 lom 


on the right. Thus to avoid ambiguities, we will use ^ to denote this variation: 


^ab _ _.abr. ,ab ‘^{df^T^f _ 12(d;.7r) 

(57r“ (fvr^ F^uj^ F'^oj'^ 

+ ^(d'^7r')d^ + d^(d^vr“)- 


(4.6) 


F2a;3 


F2^3 ■ 


Then is the correct form of the differential operator satisfying 


1 


Lr, = —TT U TT 


Aan-^ab^Ab 


Now splitting into a free part Dq^ = and an interaction part 14“^, 


jyab ^ Jjab yab^ 


we have 


yab ^^abQ 


1-io- 


[Q^. ^ jafe 2(d;.7r)2 _ 12(d^7r)^7r°7r^ 


M 9 


+ 


47r“ 




F^lo^ F^uj^ 

47r^ 


(4.7) 


(4.8) 


(4.9) 
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As in previous examples, all z-dependence is introduced by the differentials. Thus we 
decompose V°‘^{z) into two parts, 


V'^\z) - izW'^\ 

TJ/ab xab r f)— ^ ^ 1 i (a— b\ i /o— a\ 

^ =M'9 , vr ) + (5 vr 


where d = q^d^, and {•, •} is the anti-commutator. 
Furthermore at large 2 the propagator is 




52 


Az'^{d 


-i2 


^ab ^ ) 


(4.10) 


(4.11) 


From these expression, we see that since each terms of G{z) in (3.18) has fc’s (llg^(z)) ^ and 
{k + l)’s V°‘^{z), Q{z) is of order 0(z). After some careful manipulations we get following 
compact form: 


g(z) = - izw{i + + (1 + w^)-^ 


q2 

V + W 


4(5-)2 

where the second term is the boundary operator we are looking for. 




+ 0 ( 1 ) 

(4,12) 


4.2 An Alternative Decomposition of V 

In the last subsection we see that when V(z) ~ Cl(^)) the boundary operator receive 
contribution from diagrams with arbitrary number of hard propagators, which makes the 
derivation of boundary operator quite involved. In order to simplify the computation in 
this subsection we introduce a different decomposition, 


P — 2 ?^ -|- Vq, — Pq + 


(4.13) 


in which we have absorbed 0(z) terms of V to a new differential operator P^, and Vq ~ 
0(z°). 

Then G(z) can be written as^^ 


g(z) = -Vo(v-^-Vo^)Vo 
= -Po 


'0 

V 


P7^ - Pk^ Po + VoV-%(l + p-Vo)-^P,-^Po 


(4.14) 


We will assume Vz(z) ~ 0(z), then T>~^ ~ hi the second term of (4.14), PqP”^ ~ 

O(z^), and PT^Vq ~ Cl(j)- So in (1 -|- 'D~^Vq)~^ only the leading term contributes to the 
boundary operator, 

g{z) ~ - Po(P-' - Vd)Vo + PoP,-VoP-'Po ^ 

=14(1 + + (1 + ^0(1 + 

^°For simplicity in this subsection we will neglect wavefunctions ej,.. 
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To simplify this expression further, let us define 


V,iz) = V, + zX, = ^ + ^ + 

z z^ 

then 

(1 + ^ (l + rf.X + 

= (1 + diJf- i(l + diXy^d^X + rfiK,)(1 + rfiX)-‘ + 0(^-2) 

z 

After a bit of calculation we find Q[z) can be written as 


g(z) =zx(i + diX)-^ + (1 + xdi) 


\-i 


V - Xd2X 


(l + diX)-^ + 0(-) 


Now we go back to nonlinear sigma model. We will choose 


yab 

l^ab 


1 . .2 4-7r^ 47r^ 

,ab^(A^ _ 12(a;,7r)^7r°7r^ 

F2a;3 F^uj^ 


and we have 


di 



92 

4(9-)2- 


Plug (4.20) into (4.18), we find the same expression as (4.12). 


(4.16) 


(4.17) 


(4.18) 


(4.19) 


(4.20) 


4.3 A Theory with Yukawa and Quartic Scalar Coupling 

In this section we consider a massless theory with Yukawa and interactions: 

L = ^(pd'^ct) + (4.21) 

For compactness, here we will focus on the first deformation, and the discussion of the 
boundary operators of multiple step deformations can be found in Appendix A. 

The chiral fermion carries SU (2) indices, and we will use the following angle and bracket 
notations to keep track of these indices: 

I'lp) = V'a, (^1 = V''", 1^] = [^1 = (4.22) 


To find the boundary operator after the first deformation, hard helds are combined into 


= 


^ 4)\ 

lip) L = (<(> [V'l (V’l), 
ViV']/ 


(4.23) 


and the corresponding operator 11“^ can be decomposed to following free part 
tion part: 



0 

« \ 



Xi'ipl 

A(V’l\ 


0 


, = 

M^] 

X(j) 




« ) 



0 

x<p J 


and interac- 


(4.24) 
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The inverse of free part D 


OL^ 

0 


is 


{Dq )ais 


/l 0 0 \ 

0 0 1 ( 9 ^)”^. 

\0 0 ) 


(4.25) 


To analyze the z-dependence, one can simply do the replacement ^ — izq^. It is not 

hard to see that is not affected by the shift, while becomes 


{Dq )ap{z) — 


(izd~ 

0 


V 0 -F-i( 


0 

0 

a'^q-a 

2(9-)^ 




iq-u I 1_ 
a- ^ z 


0 

a'^q-a 

2(^ 

0 




+ 0 ( 4 ). 


(4.26) 


The theory has 3 physical states, therefore there are in all 9 possible ways of BCFW 
shifts. Among these shifts, {'ip\'4>] shift is 0{^), and the corresponding boundary operator 
vanishes. There are 5 0(1) shifts, and the corresponding boundary operators are 




Qim =A0(ln) + ^A2(nV’)^(nV’) 

^ CJ 


(4.27) 




+ 


^ ■ ^1 71 


using our general formula (3.18). 

The other 3 shifts are 0{z) and the corresponding G{z) operators are as follows. For 
(i?i)|'0] shift we get 


g{z) =zX[l'ip] - Xinip] + ^Xgcp'^-^ll'ip] + 


A^A / .-.q ■ a. 


'^1711 + (^l^l7lT|,fi| - 2im]^A + 0(1). 


(4.28) 


For 


shift we get 


igX 


|A|^ 


g(z) = - zX{mjj) + A(IV’) - + ^(j) ( i(ln)^[lV’] + 


, cr 




+14 


+ ("7)gr.., a 


04+2 1 ' ^-^'1 


— IV’] -2iV> ^(raV’) j +^(-) 


(4.29) 
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Finally, for (V’lV’] shift we get 

8(2) = ^ 

+ ia|2 

\ ^(y C/ 

+ ^(-^ff(nV’)^</>^^[^l] + A0(ln)^[l^-]^[^l] 

- + A[ln](nV’)^(V’?^)^</> - *A(nV>)^('0l^^I!]</>) 

+ ([^I^IV’) - (V^l^l^]) ^[^1] - 

- 2(nT/')^[V5l]^(/)^ - + 0(i) 

The boundary operators of these 3 shifts are 0(1) order of the corresponding Q(z) 
operator and they have more complicated expressions. However, in practice we can always 
choose shifts which have better large z behaviors and the corresponding boundary operators 
are much simpler. 

5 Conclusion and Discussions 

To summarize, we have introduced the boundary operator as a tool to study boundary 
contribution of BCFW recursion relations, and presented an algorithm of deriving boundary 
operators. To demonstrate our algorithm, several examples have been presented to show 
different aspects of the algorithm, like the presentation and decomposition of operator P, 
and the analysis of z-dependence. 

Another generalization of BCFW has been discussed in a recent work [18]. But instead 
of deforming external legs one by one, they choose to deform multiple legs at the same 
time. Surprisingly, these two methods have a lot of similar features. For example, any 
amplitude in 4-D renormalizable theories can be obtained by at most 5-line(or 4 step in 
our method) deformation, while amplitudes in non-renormalizable effective field theories 
usually do not have the desired large z-behavior. In the language of OPE, the multiple line 
shift corresponds to the products of more than two operators. The evaluation of boundary 
operators in this case can be much more difficult because the hard field action is no longer 
Gaussian. But still it is worthwhile to check whether it can be accomplished at least for 
simple cases like 3-line shift. 

As has been mentioned in the introduction, the boundary operator method can be used 
to extract form factors and correlation functions from scattering amplitudes. In fact, (2.5) 
tells us each 2 order of the shifted amplitudes corresponds to the form factor of a local 
operator Si. 

Last, let us point out that one of the advantages of the new description of boundary 
contributions is, since OPE can be defined at loop level or even non-perturbatively, there 
should be no obstruction to generalize boundary operators to loop level or non-perturbative 
settings. It would be interesting to study the loop corrections to boundary operators in 
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supersymmetric gauge theories, where BCFW recursion relations for planar amplitudes are 
known even at loop level [20, 21]. 
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A Boundary Contributions of Multiple Step Deformations 


In subsection 4.3 we discussed the boundary operators in a theory with Yukawa and quartic 
scalar interactions, and found only one 0{^) shift, (V’jV']- la order to compute scattering 
amplitudes with arbitrary configurations of external legs, in this appendix we study the 
boundary operator of further deformations, so we can find a simplest way to calculate the 
whole amplitude. 

Suppose the amplitude has at least 4 legs, we can use a series of {tp ■ ■ ■ 'p\'pi\ shifts to 
find the whole amplitude. From last section we find 


0(W]=A[ln](/.-^A2[l^]^[^l] 


Following the steps in subsection 3.3, we find 





= 

Hd, 

= 



V 


A[ln] 

Sfll] [;jL[,^l] + (^[,^1])] 

0 


thus we get following boundary operator 




after the second deformation, and 


(pWV'IV'l =iA 2 [^ 2 ][ 13 ]^. 


(A.l) 


(A.2) 


(A.3) 


(A.4) 


after the third deformation. Since this operator contains no soft field, it will not contribute 
to a n point amplitude if n > 4.^^ Therefore n-point amplitudes with at least 4 pj legs and 
n > 5 (or equivalently, 4 pj legs) can be computed by three recursion steps. 

If the amplitude contains at least 4 scalars, we can use a series of {4> ■ ■ ■ (j)\<p] shifts. 
After the second deformation {(p4)\4>] shift, we find boundary operator 




^ g<P 

-§iAp 

ui^i' 


= gcj) 


(A.5) 


general, if the boundary operator of a k-step shift contains no soft fields, the boundary contribution 
vanishes if n > fc. The boundary contribution does not vanish for the case n = k, unless the boundary 
operator vanishes. 
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After the third deformation we find = g, therefore the amplitude can be computed 

using three recursion steps if ra > 4 with at least four scalars. This is consistent with 
the result of the scalar theory: a pure scalar amplitude would not involve any Yukawa 
interaction, so the amplitude should be identical with amplitudes in a (p'^ theory. 

If the amplitude has at least one ip and two p legs, with a bit calculation it can be 
found that the complete amplitude can be computed using a two step recursion {cppl'ip]. 

In summary, any amplitude with n > 4 in this theory can be computed using one of 
these 4 types of shifts: {ipl'p], {(ppcplp]. Suppose 

the amplitude has least a ipip pair, then {ipltp] is a good one step shift. Pure bosonic 
amplitudes can be computed using a {(p(f)(p\(j)] shift. For other types of amplitude, without 
loss of generality, we can assume there are no p) legs. If the amplitude contains more than 4 
Ip legs, it can be computed using a {ip'tpipl'ip] shift. If the amplitude only has 2 pj legs, since 
n > 4, there are more than 2 cp legs, and the amplitude can be computed with a {'ip\pp] 
shift. 

B Yang-Mills Theory in D Dimensions 

In [19, 22], the boundary contribution in Yang-Mills theory is discussed. The analysis of 
large z behavior relies on a "spin Lorentz symmetry" of the Lagrangian aftering imposing a 
special gauge. With our new method, the large 2 ; behavior of the deformed amplitude can 
be read directly from the OPE (2.6), which can be computed explicitly. 

Let us consider the Yang-Mills theory in spacetime dimension D > A: 


r — _ 1 7?“ 

4 

(B.l) 

Splitting Yang-Mills field A^ into 


Afj, ^ Afj_ + A^ 

(B.2) 

and adding the gauge-fixing term^^. 



we find the following expression for 

= 1 ^;)“ _ 2gfahcp^.uc^^ ^Ab 

The operator F in this case is dressed with both spacetime and SU{N) indices, 

pfj.a-,ub _^iiu^abQ2 _ yi^aiub 

yl^a.ub ^ _ gJabc^2F^^yc ^ A^^}) - 


(B.3) 


(B.4) 


^^Most conventions in this subsection follow [19] and [22]. 



Denoting the inverse of the V operator as after the LSZ reduction we find 

d" = V^a;ub " + • • • (B.6) 

Now we consider the large 2 : behavior. Under the shift (l|n], we have 

= Yi^a;ub ^ 2izgQo^A^'^, (B.7) 


i.e., is of order 0{z). However, as discussed in [19] the 0(z) piece vanishes when 

light cone gauge are imposed on the soft fields^^. Most terms in (3.18) are actually at most 
O(^), except the leading term, where light cone gauge can not be imposed on the single A 
held. We will denote the leading term by Gi{z). Following the conventions for polarization 
vectors in [19] we hnd 




Gi{z) = 


yii.a-,ub^z) ^e+(fc„) e^{kn) ei{kn)J 


^Qfj, ~k zkjiij\ 

+ 2izgr^-g^^''q^A-^) {q, - zki, q, eiCkn)) 

\ / 

^Q^l “ 1 “ zkiifj\ 

-\ky + 2izgr^^g>^''q^A^^) (q, - zki, \kn. ei{kn)) 

\ / 

Gi^{z) 0 {qf, + zkn^)V^-’’'’^ei 

0 0 0 


ac\^J 


Y^V>^a-,vb^^^ _ Q Y^^y^,a■,ub ^ 2izgq^A 

where we have chosen as normalized constant vectors satisfying 


+ 0(1) 


(B.8) 


(B.9) 


From (B.8) the large z-behavior of various shifts can be read out. First (<h|g'“] and 
( 5 '*'|$] shifts with any $ are of order O(^), so boundary terms are zero. Next for {g~\g~^] 
shift, since all higher order terms in Q ^{z) are at most^^ ^(z), we hnd 

G-^{z) = g^+{z) + o{z) 

=2iz^ki • kngr^A^^q^ + ^2 [ki ■ knK^^ + 2gr'>^F^^'^^kyknu) + 0{z), 

where 

J^ab ^ ^ac| ^ ^2 jace jbde_ (g 

are imposing background gauge to high energy modes, and light cone gauge to low energy modes 
of the gauge field. Around the energy scale A, a smooth function can be used to link two different gauge 
fixing terms. 

^“^All higher order terms have at most one propagator which is O(^), while two polarization vectors 
introduce a 2 ^. 
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If we shift two transverse gluons {g^\g^] 

gig^y] =ei2izgr‘’^g>^''qaA^^ei + 0(1) = 2izg5^^ + 0(1) (B.12) 

If i = j, gyy\ = 0(z), otherwise gyy^ = 0(1). 

All the results above are consistent with those of [19]. 
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